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Abstract. We present a numerical study for the threshold percolation prob- 
ability, p c , in the bond percolation model with multiple ranges, in the square 
lattice. A recent Theorem demonstrated by de Lima et al. [B. N. B. de Lima, 
R. P. Sanchis, R. W. C. Silva, STOCHASTIC PROC APPL 121, 2043-2048 
(2011)] states that the limit value of p c when the long ranges go to infinity 
converges to the bond percolation threshold in the hypercubic lattice, 7L d , for 
some appropriate dimension d. We present the first numerical estimations for 
the percolation threshold considering two-range and three-range versions of 
the model. Applying a finite size analysis to the simulation data, we sketch 
the dependence of p c in function of the range of the largest bond. We shown 
that, for the two-range model, the percolation threshold is a non decreasing 
function, as conjectured in the cited work, and converges to the predicted 
value. However, the results to the three-range case exhibit a surprising non- 
monotonic behavior for specific combinations of the long range lengths, and 
the convergence to the predicted value is less evident, raising new questionings 
on this fascinating problem. 

1. Introduction 

Percolation is one of the most studied phenomena in statistical physics since it is 
closely related to several phase transitions in a wide range of systems. As a typical 
interdisciplinary subject, the literature shows contributions of researchers from dif- 
ferent areas as mathematics, physics, civil engineering, soil science, hydrology etc. 
treating a wide range of problems from infiltration and porous flow to jamming 
and glassy transitions, passing by order-disorder models and magnetism [TJ [2]. A 
definition to percolation could be simply stated as following: considers a regular 
two-dimensional square lattice (Z x Z ) in each the vertex are connected by bonds 
which are present with probability p; with we start with one single bond occupied 
and randomly sort a new bond to connect two neighbors sites, eventually a large 
cluster will span along the entire system and we say that the system has percolate. 

Despite its simplicity, several open problems still challenge any theoretical treat- 
ment; for example, the exact value for the threshold percolation probability, p c , 
for the regular lattice in three dimensions is not known, and even the continuity 
of the probability density function at p c still is an open problem ( for general facts 
in percolation see [3J. The major importance in determine the percolation thresh- 
old and how it depends on the dimensionality of the space, D, is due to the link 
between percolation and phase transitions in natural systems. Several emergent 
phenomena arise near to the critical point, where power-law and fractal behaviors 
are observed. Most of these phenomena could be explained considering the cluster 
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size distribution of similar components along the system, and percolation plays a 
major role to understand the mechanism involved in these natural processes. 

In the recent years, a particular attention was paid in literature on the study of 
irregular lattices, or complex networks. A central characteristic of these networks 
is an inhomogeneous distribution of bonds, or links, by node, and even a wide dis- 
tribution of bond lengths. Hence, a natural question arises: what is the percolation 
threshold and how it behaves in function of the network parameters? Several works 
[5] have studied this problem, but most of them have focused the applications, e. g., 
the Ising model in complex networks [6 . Here, we present a more generic study of 
this problem from the percolation theory point of view; we introduce a percolation 
model with multiple bond ranges extending the usual bond percolation definition 
in the square lattice: we draw a bond to connect two sites in the lattice randomly 
chosen among a list of available sites within a given range. Thus, by means of this 
generalization of the standard model we expect to unveil possible hidden features 
which arise when we combine several bond ranges, and expect to shed some light 
in the problem of percolation in more complex networks. 

Our primary aim is to verify a very recent conjecture from de Lima et al [3] which 
states that the percolation threshold should be a decreasing function of the second 
range, when two bond ranges, (l,fc), are available, and Theorem 1 in [3] which 
states that the p c value in this case converges to the bond percolation threshold in 
a d = 4 hypercubic lattice in the limit of k going to infinite. It worth to mention that 
add a new bond range is equivalent to increase the freedom degrees of the system, 
thus an operation analogous to increase the spatial dimension D. Here, we present 
an extensive numerical study of the percolation threshold, p Cl for two-range, (1, k), 
and three-range, (1, k, m), bond percolation models. In the three-range model, the 
predicted value when k — > oo and (m/k) — ¥ oo, is the percolation threshold of the 
d — 6 hypercubic lattice. Surprisingly, the numerical results have shown a non- 
monotonic behavior for specific combinations of k and m which are not predicted 
by the theory. 

The paper is structured as following: in order to well state the problem and 
present the analytical conjectures, we present the model definition and the theoret- 
ical predictions in the limit of diverging ranges in the next section. The numerical 
procedure and simulation results are discussed in the following and we drawn some 
conclusions and perspectives at the final section. 



2. Theoretical Aspects 

Let G = (V, E) be a graph with a countably infinite vertex set V. Consider the 
Bernoulli site percolation model on G which associates to each vertex the values 1 
( "occupied" ) and ( "vacant" ) with probability p and 1 — p respectively. This can 
be done considering the probability space (fi, F, P p ), where f2 = {0, 1} V , F is the 
er-algebra generated by the cylinder sets in f2 and F p = Yl v<£¥ ^( v ) * s the product 
of Bernoulli measures with parameter p. We denote a typical element of fl by u. 

Given two vertices v and u, we say that v and u are connected in the configuration 
w if there exists a finite path (v = vq, v±, . . . , v n = u) of occupied vertices in V, such 
that Vi Vj, Vi ^ j and (vi,Vi+\) belongs to E for alH = 0, 1, . . . , n — 1. We will 
use the short notation {v 4-> u} to denote the set of configurations where u and v 
are connected. 
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For a given vertex v, the cluster of v in the configuration uj is the set C v (ui) — 
{u G Y;v -H- it on w}. We say that the vertex v percolates when the cardinality 
of C v (u>) is infinite; we will use the following standard notation {v <H- 00} = {u; e 
f2; #C v (ui) = 00}. Once v fixed, we define the percolation probability of the vertex 
v as the function 8 v (p) : [0, 1] 1— ► [0, 1] with 9 v (p) = P p (v O 00). The percolation 
threshold (or critical point), p c (G), is defined by 

p c (G) = swp{p€ [0,1];%) = 0} . 

From now on, the vertex set V will be Z d , d > 2 and for each positive integer k 
define 

E k = {((«!, ... ,u d )) eVxV;3!ie{l,...,d} 

such that \vi — Ui \ — k and Uj = v>j, V j 7^ i} . 

Let's define the graph G fe = (V, Ei UE^) that is, G k is the 1 d equipped with nearest 
neighbors bonds and long range bonds with length k parallel to some coordinate 
axis. Observe that G k is a transitive graph, hence the function v (p) k does not 
depend on v and we write only 9 k (p) to denote P„(0 -H- 00) on the graph G k . 

Consider the sequence, (p c (G k ))k, of the percolation thresholds of the graphs 
G k . In [3] (see Theorem 1), it is proven that 

lim p c (G k ) -p c (Z 2d ),Vd> 2 , 

k— ^+00 

that is, the percolation threshold of G k converges to the percolation threshold of 
the hypercubic (nearest neighbor) lattice in 2d dimensions. It is also conjectured 
[3] that the sequence {p c {G k ))k is non increasing in k. One goal of this work is 
to simulate long range percolation on the graph Gk in d = 2. These simulations, 
as described in the next section, show that p c (G h ) goes to p c (Z 4 ) monotonically, 
confirming the Theorem 1 and the Conjecture in 

Theorem 1 and the Conjecture in [3] can be generalized for bond percolation 
with multiple ranges. Given a sequence k = {k\, . . . , k n ) with ki £ {2, 3, . . . }, Vi, 
consider now the graph G^ as (Z d ,Ei U (U™ =1 E fclX ... xfe J). That is, G~^ is Z d 
decorated with all bonds parallel to each coordinate axis with lengths l,hx,ki x 
fc 2 , . . • , fci x fc 2 x ••• x k n . In this case Theorem 1 of [3] states that 

lim Pc (G k ) = Pc (Z d{n+1) ),W > 2 

k-i-^ooyi 

and the Conjecture says that the probability of percolation on the graph G k is non 
decreasing in each variable ki. 

In this work, we also perform simulations in the case n = 3, considering bonds 
with three different ranges 1, k and m. The results support the Theorem 1 and 
the Conjecture in [3] (in this case, m should be a multiple of k). One interesting 
feature displayed by simulations is that the percolation threshold is not monotone 
in k, with m fixed. 

3. Numerical Procedure 

In this section we present the numerical approach used in the simulations. We 
consider a square two-dimensional regular lattice as the site substrate of the model. 
When a link to a given range is open, two sites in a column, or in a line, will be 
connected to each other forming a bond - see Figure [T] For each bond, besides to 
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Figure 1 . Neighborhood of a bond in the 2-range version of the 
model. The first neighbors of the central black bond are shown, for 
(l,fc) = (1,2). Note that there is 6 short-range bonds (A-F) with 
the same length as the central bond, and 8 long-range bonds (1-8). 
The sites connected to the central bond are filled with different 
gray levels, depending on the range of the bond linked to the site. 



the six neighbors in the same range, there are eight neighbors consisting of bonds 
of another length. Thus, in the 2-range version of the model there are 14 neighbors 
and, in the 3-range, 22 neighbors to be taken into account for every new bond 
opened. 

The algorithm used to estimate the percolation threshold is based on the work 
of Newman and Ziff [7J . Basically, we label each bond of the lattice assigning the 
first 2N natural numbers to bonds of unitary length (where N = L x L, and L is 
the system size); the labels 2N + 1 to 4N are assigned to the range k, and so on. 
Thus, we enumerate all the possible bonds to be opened in the system from to 
r 2N, where r is the number of ranges considered. 

In this way, each realization consists of a random permutation in the list of the 
bonds to be opened, varying the occupation probability, p, from to 1. After 
each bond is open, the percolation condition is tested using the modified Hoshen- 
Kopelman algorithm to identify the clusters pp. The probability of a given site to 
belong to the infinite cluster is then estimated from the simulated data using the 
normalization procedure described in [7J. 

The code was developed in C and Purebasic ® languages. Purebasic ® code 
developed by M. Schnabel is procedure-oriented and 1000 samples of a system with 
L = 1024, single range percolation model, takes around two minutes to run in an 
AMD Athlon X64 1.9 GHz and 3 GBytes DDR2 RAM machine. It is about 20% 
faster then the similar C implementation. 
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Figure 2. Probability of a site belongs to the infinite cluster in a 
single range bond percolation model. The plot shows the variation 
of to 6 values of L = 32, 64, 128, 256, 512, 1024. Note that all 
curves tend to cross at the same point. Considering two consec- 
utive system sizes, we measure the ordinate value of the crossing 
point. The inset shows the crossing values obtained in function 
of the inverse of the geometric mean between the system sizes, 
1/ y/LiLj . The independent term of the linear fit correspond to 
the best estimative for the percolation threshold to a infinite sys- 
tem size, and matches perfectly with the exact value p c = 1/2. 



In order to validate the code developed, we initially consider the bond percolation 
model limited to a single range, in a square lattice of linear extension L, with 
periodic boundary conditions in both directions. We simulate 1000 realizations 
for each system size, 16 < L < 4096, and calculate the probability of a given site 
belongs to the infinite cluster, Figure[2] In Figure[3]we show the corresponding 
calculation for the two-range bond model considering k = 7. We denote by (l,k) 
or (l,k,m) the 2-bond and 3-bond models, respectively, where k, m G N. 

The percolation criterion used assumes that the infinite cluster arises when there 
is a path of connected bonds crossing the entire system, vertically and horizontally, 
not using the periodic boundaries. Although the particular choice used for this 
criterion, we believe that we would get qualitatively identical results even for other 
choices, as discussed in the next section 
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Figure 3. Probability of a site to belong to the infinite clus- 
ter in a 2-range (l,fc) bond percolation model, with k = 7. 
The plot shows the value of $(p) considering 6 values of L — 
32, 64, 128, 256, 512, 1024. The inset shows the finite size analy- 
sis used to obtain the best estimative for the percolation threshold 
at the infinite system (see the caption of Figure [2| . This proce- 
dure is repeated for each value of k considered in this version of the 
model, and the dependence of the threshold percolation in function 
of k is shown in Figure [4] . 

4. Results and Discussion 

Figure [3] presents the procedure used estimate the percolation probability to the 
2-range model (l,fc), with k — 7. Considering different values of A;, it is straight- 
forward the calculation of p k c for increasing k, which is shown in Figure |Zj The 
limiting value for k — > oo is p c ~ 0.162(2), close to the expected value for the 
bond percolation threshold at the d = 4 hypercubic lattice - p c = 0.160131(1) [5] as 
stated by Theorem 1 [3J. Thus, we can conclude from this plot that the percolation 
threshold exhibits a monotonic dependence on k, supporting the Conjecture in [3J, 
and confirming the prediction for the limit k — > oo stated by Theorem 1. 

The results for the 3-range model are shown in Figure [5] We consider five 
different values to m, and all the possible k between 1 and m. All results are shown 
in Figure [5] We can observe some striking features: first, there is a clear downward 
trend of p c with increasing to. To to fixed, generally p c decreases with k, except 
when k approaches to m, or to the half-value of to. This trend was observed for 
to even, odd or prime. For to = 32 a curious feature appeared: the value of p c 
increases about k = 11 and k = 21. We do not have a good explanation for this 
match so far. The increase to k = 8 can be also be attributed to the correspondence 
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Figure 4. Threshold percolation probability, p c , for 2-range bond 
models, with increasing range k. There is no justification for the 
fit function used and the line is only to guide the eyes. The limit 
value of p c is obtained from the fitting parameter yO. 



with m = 32, as 32 is multiple of 8; however, for ranges of k — 4 and to = 16, for 
example, a similar effect was not observed, discarding this conjecture. 

Considering the curves for m = 10, 13 and 16, they all behave similarly, with p c 
decreasing with k except for k around to to or to/2. For to = 13, the increasing at 
to/2 is more discreet, but occurs between k = 6 and 7, since the half- value of to is 
not an integer. The curves for to = 32 and 33 show different features. We observe 
that the decrease of p c with k is more subtle, but the increasing at to and to/2 
is quite evident. The peaks observed at k = 8, 11, 21 are unexpected and not well 
understood for the moment. The limit value of p c ~ 0.0955(7) around to k = 25, 
approaches to the bond threshold percolation value at d = 6 hypercubic lattice, 
p c = 0.0942019(6) [8], which is expected when the two long ranges tend to infinity. 

5. Conclusions and Perspectives 

We present by the first time the percolation threshold estimation for a multiple 
range bond percolation model. We study the behavior of p c in two and three 
range models in function of the bond range using numerical simulations. We have 
confirmed the predicted value to the percolation threshold when the long ranges 
tend to infinity in these models - Theorem 1 in [3] . A remarkable effect was observed 
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Figure 5. Threshold percolation probability, p c , for 3-range bond 
models, in function of the range k; five different curves correspond- 
ing to increasing m are shown. Note the differences when m is 
odd or a prime number, smoothing the curves. The striking non 
monotonic behavior observed for specific combinations of the long- 
ranges still is not fully understood to the moment. The dashed 
line represents a Lorentzian fit to the minima of the plots. 

for the 3-range version of the model, with an unexpected increase in the percolation 
threshold for some values specific values of the second and third ranges which defy 
any convincing explanation to the present date and could represent an important 
feature in models with long range interactions or even systems dealing to complex 
networks. 
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